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ABSTRACT 

This article presents the first direct “observation” of the global-scale, 3D 
coronal magnetic fields of Carrington Rotation (CR) Cycle 2112 using vector 
tomographic inversion techniques. The Vector tomographic inversion uses obser¬ 
vational measurements of the Fe XIII 10747 A Hanle effect polarization signals 
by the Coronal Multichannel Polarimeter (CoMP) and coronal density and tem¬ 
perature structures derived from scalar tomographic inversion of STEREO/EUVI 
coronal emission lines (CELs) intensity images as inputs to derive a coronal mag¬ 
netic field model that best reproduces the observed polarization signals. While 
independent verifications of the vector tomography results cannot be performed, 
we compared the tomography inverted coronal magnetic fields with those con¬ 
structed by magnetohydrodynamic (MHD) simulations based on observed photo- 
spheric magnetic fields of CR 2112 and 2113. We found that the MHD model for 
CR 2112 is qualitatively consistent with the tomography inverted result for most 
of the reconstruction domain except for several regions. Particularly, for one of 
the most noticeable regions, we found that the MHD simulation for CR 2113 
predicted a model that more closely resembles the vector tomography inverted 
magnetic fields. In another case, our tomographic reconstruction predicted an 
open magnetic field at a region where a coronal hole can be seen directly from a 
STEREO-B/EUVI image. We discuss the utilities and limitations of the tomo¬ 
graphic inversion technique, and present ideas for future developments. 
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Introduction 


Routine direct measurements of the coronal magnetic fields are indispensable for our 
eventual understanding of the nature of solar coronal phenomena at all spatial and temporal 
scales. I mportant progresses were made in the last two decades in t he development of the 


theories (jCasini & Juc 


observing techniques flLin et al. 


m 


1999 


2000 


in fc Casinill2000l:ICasini fe Lirul20021) . instrumentation and 


2004J) to directly measure the polarization of coronal 


emission lines (CELs) that bear the information of the coronal magnetic fields. A major mile 
stone was reache d w ith the deployment of the Coronal Multichannel Polarimeter (CoMP, 
Tomczvk et al.l 120071 . 120081 ). CoMP provides synoptic linear polarization measurements of 
the Fe xill 10747 A emissions which directly maps the orientation of the magnetic fields 
of the solar corona. However, due to the low density and opacity of the solar atmosphere 
the observed polarization signals are the total signals summed along the line-of-sight (LOS) 
of the measurements with non-uniform temperature, density and magnetic field structures. 
Therefore, except in a few special cases, direct inference of the 3D coronal magnetic field 
structure from polarization data is in general not possible. 

Parallel to the development of the observing capabilities, a vector tomography inversion 
method to infer the coron al magn etic fields based on CEL circular polarization data was 
developed by iKramar et al.l 020061 ). More recently a second vector tomography inversion 
method that combines coronal linear polarization data and information of the temperature 
and densit y structure o f th e co rona to reconstruct the 3D coronal magnetic fields was de¬ 
veloped by IKramar et ah) (120131 ). We applied this technique to CoMP data obtained during 
Carrington Rotation (CR) cycle 2112 to derive the global 3D coronal magnetic fields up to 
1.3 -Rq. However, as it was pointed out in Kramar et ah (12013)) . while polarizations of the 
CELs respond directly to the local magnetic fields, their visibility, i.e., the strength of the 
emission lines, strongly depends on local thermodynamic properties of the plasma.Therefore, 
3D temperature and density structures of the corona are also needed as inputs for the vector 
tomographic inversion program for the calculation of the strength of the coronal emission 
lines. In this work, the 3D temperature and density structure of the corona were derived by 
scalar tomographic inversion of the ultraviolet coronal images obtained by Extreme Ultra¬ 
violet Imager (EUVI) on board of the Solar Terrestrial Relations Observatory (STEREO) 
spacecraft. In the following sections, we will first describe the scalar inversion we used to 
derive the temperature and density distribution of the corona from STEREO/EUVI data. 
Then we will present results from the vector tomographic inversion of the 3D magnetic fields 
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of CR 2112. Limitations of this methods, as well as future improvement of the technique 
will be discussed. 


2. Reconstruction of Density and Temperature 


The process of the 3D reconstruction of the coronal electron density and temperature 
based on STEREO/EUVI data consists of two steps. First, the 3D reconstruction of EUVI 
emissivities. Second, finding the density and temperature values which satisfies the pre¬ 
viously reconstructed emissivities for every voxel in the 3D reconstruction domain. While 
first step is t he same as in the Differen t ial Emission Measure Tomography (DEMT) method 
described bv lFrazin et al.l (120051 l2009h : IVasquez et al.l feOlOl h our implementation of second 
step is different. The main differences in second step are that, in our method, we assumed 
the plasma is isothermal within each volume element (voxel) and different implementation 
of an inversion procedure. 

Briefly, the measured coronal emission I\ in any of the EUVI channels can be represented 
by the LOS integral 

h = k- f £\(r)d£, ( 1 ) 

LOS 

where £\{f) is the emissivity at the position r in the corona for the selected channel A, i.e., 
light intensity (in photons per second) emitted per unit volume, per unit solid angle. The 
coefficient k accounts for pixel size, aperture, distance to the Sun. The inversion _of [1] for 


£i 7 i (r ), £195 (r), and £ 284 (u) is performed by tomography method described in iKramar e t a h 

ihoiA 


For an locally isothermal plasma, £\(r) can be expressed by (Phillips et al. 20081 ) 


£ X = G(T,N e )Nl 


( 2 ) 


where T is the temperature, N e is electron density (and a function of T), and G(T,N e ) 

For our inversion G(T , N e ) was calculated with the CH1- 


is the contri 
ANTI code 


oution func 


Dere et al 


ion. 


20091 ). and the Fe ions abundances were derived from results of 


Arnaud fe Raymond! f 1992^ 7 For every grid cell, we derive three (one for each EUVI chan¬ 
nel) N e {T ) curves from ([2j) based on the emissivities £\{r) derived by the tomography from 
equation ([T]). With the local isothermal approximation, the three N e (T ) curves will converge 
at a point indicating the solution for T and N e for a given grid cell. 

We applied this inversion method to the STEREO/EUVI ( Wuelser et al.ll2004 ; Howard et al 


20081 ) data of CR 2112. STEREO/EUVI observes the coronal up to approximately 1.7 Rq 






























4 


in four spectral channels (171, 195, 284, and 304 A). We used the off-limb data of images 
from the 171, 195, and 284 channels for the tomographic inversion of the temperature and 
density. Similarly to Vasq uez et ah (120101 ). three average images per channel per day over 
half a solar rotation during CR2112 were used for the reconstruction. The average images 
were obtained by averaging three off-limb images taken with a two hours interval. The res¬ 
olution of the EUVI images were also reduced to 512 x 512 pixels. The upper boundary 
of the reconstruction domain is limited to 1.57?© due to the lower signal-to-noise ratio of 
data higher in the corona. The lower boundary i s set to 1.057?^ in o rder to avoid opti¬ 
cal depth effects near the limb (IFrazin et al.l 120091: lYasquez et al.l 120101) . Figure [T] shows a 
spherical cross-section of the 3D temperature and electron density for CR 2112 at 1.177?© 
reconstructed with this method. Additionally, the density map at 2 .07?^ obtained by tomo¬ 
graphic inversion of the STEREO/COR1 data (Kr arnar et al. 2009) are shown. This map 
will be used in later sections to verify the magnetic field of the corona. 


3. Reconstruction of Magnetic Fields 


The CoMP instrument measures the intensity and the linear and circular polarizations 
(Stokes I,Q,U,V) of the Fe XIII lines at 10747 and 10798 A of the solar corona from ~ 1.03 
to ~ 1.47?©. In the present research, we use data only from the 10747 A forbidden line. 
Averaged daily CoMP Stokes Q and U images from July 13 to July 26, 2011, corresponded 
to CR 2112, are used for this study . The re constr uction of the coronal magnetic fields 


follows the method described in [Krarnar et ah (1201311 . using the 3D T and N e derived from 
tomographic inversion of STEREO/EUVI data as fixed input parameters. Although T and 
N e were derived out to 1.57?©, examination of the results indicated that the inversion above 
1.37?© may not be very reliable. Therefore the outer boundary of the reconstruction domain 
for the magnetic fields were limited to 1.37?© for this work. 

Figure [2] shows the LOS projections of magnetic hied lines originated at Carrington 
longitudes ±90° from LOS directions of 110.5, 120.5, and 200.5° during CR 2112. In this 
figure, green lines represent reconstructed held lines derived from the vector tomographic 
i nversion based on CoMP data, and the red dashed lines correspond to a MHD model 
( Mikic et al. 2007 ; Lionello et al. 2009 ) NR that was used as the starting held for the inversion. 


Grayscale areas in the background within [1.06; 1.28]7?© show the reconstructed electron 
density. Maps of the whole reconstructed 3D helds are shown in the on-line version of this 


Rww.predsci.com 

2 The MHD model used as our initial field may not be the latest result on the PredSci website. 
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paper. 

The jagged appearances of some of the field lines in the reconstructed model can be 
attributed to errors in 3D density and temperature reconstructions, the decreasing magnetic 
field strength as a function of height which makes the regularization less effective at larger 
heliocentric distances, and reconstruction errors near outer boundary of the reconstruction 
domain due to limited field of view (FOV, 1.3_R Q ) of the input polarization data. Neverthe¬ 
less, the reconstructed fields possess sufficient details to analyze the large scale structure of 
the coronal magnetic fields. It is interesting to note that the reconstructed fields are in gen¬ 
eral qualitatively similar to the starting fields except in some regions such as those labelled 
as A, B, and C, in panel 2(a), 2(b), and 2(c) respectively, where the tomographic inversion 
resulted in an open field while the starting fields were in a closed configuration. 


While the visual comparison of the magnetic field configurations between the MHD and 
tomographic inversion models in Figure [2] are instructive, a more quantitative presentation 
of the model differences is shown in Figure [3] and Table [0 In the figure, the histograms 
of the distributions of the angles, 5 , between the magnetic field vectors derived from the 
two methods in volumes within [1.05; 1.25] and [1.10; 1.25]i? 0 in heights, and within ±10° in 
longitude and latitude around regions A, B, C, and D are plotted in conjunction with that of 
the entire solid angle volume of 47T (curve E). Region D, as shown in Figure [2][b), is a region 
in which the two models agree fairly well. Therefore, its angular difference histogram serves 
as a measure of the discrepancy between the MHD and tomography models. We found that 
the distribution of the angular difference in the reference region D peaks around 14°. On 
the other hand, the magnetic field orientation difference between the two models are larger 
in regions A, B, and C, with their histograms having several peaks above 20°. Histogram 
E have two peaks at about 0 and 20°. The peak around 0° (no change in the field) is 
predominantly due to contributions from voxels below l.li? Q . This could be because the 
real field indeed very close to the starting field at lower heights, or the tomography method 
is less sensitive close to the inner boundary of the reconstru ction domain a s obse rved earlier 
in scalar field tomography for the coronal electron density (IKra ma r et ah 120091 ). Standard 


deviations, a = \J/N, of S for the regions of interest are summarized in Table [U 
Region D has the smallest value of a demonstrating fairly well local ag reement betwee n the 


two models. Region C may contain small stream er or pseudo-streamer flZhao fe Webbl 12003 


Wang et al.ll2007l: IWang||2015l: lAbbo et al.ll2015l) where the sign of the field radial component, 


B r , changes abruptly in the angular dimensions causing the large value of a. Field differences 
in region A are predominantly confined within narrow longitudinal range compared to the 
region width. Indeed, cross-sections for 10.5 and 30.5° exhibit a fairly well local agreement 
between the two models near region A. This explains a comparable values of a for regions 
A and E. Region B exhibits the same sign of B r causing lower value of a. 
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Although we have used the MHD simulated model as our reference model for evaluation 
of the efficacy of the tomographic inversion method, at this point it should be stressed that 
it is currently impossible to perform independent verification of the tomographic inversion 
results — MHD models based on observed photospheric magnetic fields inputs were derived 
with a completely different set of assumptions, and are subjects to its own limitations and 
uncertainties, and at this point which method provides more accurate result cannot be 
assessed. However, there are other coronal observations and reconstructions that in certain 
situation can help us assess the validity of the tomographic inversion results. For examples, 
coronal holes should be associated with open field regions, and current s heet structures ar e 
known to be associated with closed dipole-like magnetic field structures (iK ramar et al. 20141 ). 
Therefore we have examined the 3D electron density structures of CR 2112 at 1.17i? 0 and 
2.0-Rq, as shown in Figures[T](b),[21 and[T](c), to see if they can provide additional information 
to help us assess the validity of the tomographic inversion results. 


For region A, the density map at 2 R & showed a small low density region, inconsistent 
with the existence of a current sheet with a underlying closed field region. Thus, the open field 
configuration of the tomographic inversion still cannot be ruled out. Also, it should be noted 
that the photospheric magnetic fields used as the boundary condition of the MHD model 
contain data that were collected about 3/4 solar rotation before the reference date of the 
region A of CoMP observation, and may represent an older field configuration. Examination 
of the MHD model for the next CR 2113 as shown in Figure [2](d) reveals that it has a 
magnetic field configuration similar to that of the tomography inverted results in region A. 


For region B, the density map of this region at ~ 1.2 R Q shows a small low density 
region nearby, while the density map at 2.0 R Q shows an extended low density region (Figure 
[Q, strongly suggesting the presence of a coronal hole at this location. Additional support 
for this interpretation can be found in the EUVI 195 A images of STEREO B. During July 
2011, the angle between STEREO B and the Earth was approximately 90°. Therefore, the 
east limb in CoMP observations are observed simultaneously in the central meridian in the 
STEREO EUVI images. Figure [4] shows EUVI-B 195 A image with a central meridian 
located at 4>los = 33°. It shows a coronal hole at longitude of 30° that spread up to latitude 
of at least 50°, providing further support for the open field configuration derived from the 
tomographic inversion. 

For region C, the density cross-section in Figure [0(b) shows a narrow structure run- 



scale streamer cores may be comparable to the resolution of our tomographic reconstruction, 
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Figure [2jc) indeed shows a small closed loop near the limb embedded inside an open field 
structure, suggestive of the existence of a small scale streamer or pseudo-streamer. 


4. Discussions and Conclusions 


We have performed the first “direct observation” of the 3D coronal magnetic field using 
a combination of vector tomographic inversion of IR CEL linear polarization data and scalar 
tomographic inversion of CEL intensity data in the UV wavelengths. We compared the 
tomography inverted coronal magnetic fields with those constructed by MHD simulation 
based on observed photospheric magnetic fields of CR 2112 and 2113 as well as with the 
STEREO/EUVI 195 A image and with the global 3D coronal electron density structure 
obtained by tomography based on STEREO/COR1 data. This paper present a first look 
and preliminary analysis of the results. A comprehensive analysis of the 3D coronal fields 
and its relationship to the temperature and density structures, as well as comparison with 
other coronal models will be presented in a future paper. 


Our analysis so far suggest that tomographic inversion based on CEL polarization data 
and MHD simulations based on photospheric boundary conditions are both capable of de¬ 
riving the large-scale coronal magnetic field structures, albeit they are subject to different 
constraints and uncertainties and possess different type of errors. For examples, global 
extrapolation methods such as global MHD, potential field (PFSS), and force-free field ex¬ 
trapolations use photospheric fiel d data co lle cted over a solar rotation period as a boundary 


condition (see Section 1 in K r arn ar et ah (120141 1 and references therein) while the vector 


tomography primarily uses direct coronal data that requires observations from only half of 
a solar rotation period. Furthermore, in the tomography, the coronal fields are retrieved 
without any presumption (such as the potential or force-free conditions) on the conditions of 
the fields. The tomography requires simultaneous observations of the same object through 
multiple sight lines. Nevertheless, although current instrumentation capabilities do not allow 
for simultaneous observation of the CEL polarization from more than one sight line, data 
obtained over an extended period of time as the Sun rotates can be interpreted as data from 
multi-sight-line observations for slowly evolving coronal magnetic field structures. Therefore, 
the coronal magnetic field structures derived from the tomographic inversion method should 
be a fair representation of slowly evolving, large scale magnetic structures. 


Depending on the positions of a coronal structure relative to an observer, many regions 
close to equator are visible during a period less than half a solar rotation. Therefore, al¬ 
though the tomography method is limited in reconstructing a Coronal Mass Ejection (CME) 
structure, the method can be used for reconstruction of pre- and post-CME global 3D coro- 





nal structures flKramar et al.l l201ll ) . This will help us to understand properties of CMEs and 
their ultimate origin and precursors. 


This research utilizes polarization data from only one spectral line. To further improve 
the accuracy and reliability of this technique, we plan to incorporate additional polariza¬ 
tion data from the Fe XIV 5303^A^ line from the visible CoMP instrument at Lomnicky 
Stit Observatory in Slovakia ( Rvbak et al. 20101. We also plan to include STEREO/EUVI 
on-disk data to improve the accuracy of the temperature and density inversion. Finally, the 
current vector tomographic inversion utilizes only CEL linear polarization data which only 
yield information about the direction of the magnetic field projected in the POS. As the 
sensitivity of the instruments continue to improve, incorporation of the CEL line intensity 
and circular polarization data that carry information about the height of the data sources, 
and the strength of the longitudinal component of the magnetic fields will further improve 
the accuracy and sensitivity of the vector tomographic inversion technique. Future projects, 
such as the Daniel K. Inou ye So lar Telescope currently under construction on the summit 
of Haleakala flElmore et al.l 12 0141 ). or the Coronal Solar Magnetism Observatory (COSMO) 
proposed by the High Altitude Observatory and its partners, will be able to provide the 
full-Stokes data for the full-vector tomographic inversion, leading to a more accurate 3D 
reconstruction of the coronal magnetic field. We envision that with improved observing 
capabilities in the future, in particular, with the implementation of multiple sight line ob¬ 
serving capabilities similar to the STEREO mission, tomographic inversion method should 
provide the most reliable coronal magnetic field data to validate other coronal field inference 
methods and theoretical models. 
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Fig. 1.— The spherical cross-section of the reconstructed 3D temperature (top) and electron 
density (middle) at heliocentric distance of 1.17 -Rq obtained with EUVI Tomography, and 
electron density map at 2.Oi? 0 obtained with STEREO/COR1 tomography (bottom). 
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Fig. 2.— Panel (a), (b), and (c) show the magnetic field lines of the coronal fields of 
CR 2112 derived from tomographic inversion (green lines) and the MHD model (red lines) 
used as the initial field configuration of the tomographic inversion, originating from Carring¬ 
ton longitudes ±90° of 110.5, 120.5, and 200.5, respectively. The gray filled areas within 
[1.06; 1.28]i? 0 show the reconstructed electron density structure in these cross-sections. Blue 
lines mark boundaries between regions with open and closed magnetic field lines for the 
MHD model. The area between the inner and outer dashed circles indicate a FOV region 
where CoMP input data are used. The animation with a full set of Carrington longitudes is 
available in the electronic supplemental material accompanying this paper. Panel (d) shows 
magnetic field lines of the MHD model for CR 2113 originating at Carrington longitude of 
20.5° and 200.5°, as in panel (a). 
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Fig. 3.— Distributions of the angular differences between the tomography reconstructed 
magnetic field vectors and that of the MHD model used as the starting held of the tomo¬ 
graphic inversion in regions A (green), B (dashed red), C (black), D (black dashed), and E 
(blue). The distributions are within ±10° in longitude and latitude from longitude and lati¬ 
tude of (20.5°;-40.5°) (region A), (30.5°; 46.5°) (B), (290.5°; 40.5°) (C), and (30.5°;-20.5°) 
(D). Region E covers entire solid angle of 47T. The distributions are within heliocentric 
distances of [1.05; 1.25] (left) and [1.10; 1.25]i? Q (right). 



Fig. 4.— STEREO-B/EUVI image taken on July 20, 2011 when CoMP POS was at the 
longitude of 33°. Central meridian of the EUVI image is located at 33°. White dots show 
Carrington coordinate grid for every 30°. 
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Region 

A 

B 

C 

D 

E 

[1.05; 1.25]R 0 

a 

36.4 

32.0 

43.5 

23.2 

38.3 


24.0 

30.0 

43.0 

14.0 

0.0 

[1.10; 1.25]i? 0 

a 

39.4 

35.7 

45.4 

24.2 

41.4 

8 m 

24.0 

30.0 

43.0 

14.0 

19.0 


Table 1: Standard deviation, a, and position of the histogram maximum value, S m , for the 
analyzed regions. 
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